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Abstract 

Till date, the only consistent description of the deconfinement phase of the Sakai-Sugimoto model 
appears to be provided by the analysis of [Tj. The current version of the analysis, however, has a subtlety 
regarding the monodromy of quarks around the Euclidean time circle. In this note, we revisit and resolve 
this issue by considering the effect of an imaginary baryon chemical potential on quark monodromies. 
With this ingredient, the proposal of [Tj for investigating finite temperature QCD using holography is 
firmly established. Additionally, our technique allows a holographic computation of the free energy as a 
function of the imaginary chemical potential in the deconfinement phase; we show that our result agrees 
with the corresponding formula obtained from perturbative QCD, namely the Roberge-Weiss potential. 


1 Introduction 


The Witten-Sakai-Sugimoto model mm has been of profound significance in investigating QCD using 
holography jit]. Witten’s model for a holographic dual of pure Yang-Mills theory [2] involved essentially N c 
D4 branes on R 3 x a Scherk-Schwarz circle. Sakai and Sugimoto [3j added external fundamental quarks 
to this model by adding Nf D8 and D8-branes as probes and demonstrated that this model explains the 
low energy dynamics of QCD such as chiral symmetry breaking and the meson spectrum. By now, a large 
body of literature has grown (see, e.g, 00 for recent developments) which successfully explains various 
qualitative aspects of QCD at strong coupling using this model. 

It is, thus, natural to apply this model to systems at finite temperature. Such a study, indeed, was 
initiated by Aharony et al [7], who, following |5], employed the black D4-brane geometry as the gravity 
dual to the deconfinement phase in QCD and explained chiral symmetry restoration at high temperature. 
This was followed by many authors who investigated this model to explore various outstanding problems 
in thermal QCD such as the phase structure including the finite baryon chemical potential, which is 
difficult to compute in lattice gauge theory due to the well-known sign problem ® maim nans]. 

Ref. J7], however, also pointed out that the correspondence between the black D4-brane geometry 
and the deconfinement phase in QCD is not without problems, since there must be at least one phase 
boundary between these two phases, which, therefore, are not smoothly connected. This issue was further 
explored by two of the present authors (MM) [T], who showed that the black D4-brane geometry cannot 
describe four dimensional QCD due to issues related to the large-A), volume independence mm (as 
applied to the Scherk-Schwarz circle). Instead, [1] proposed that the deconfinement phase in QCD is 
represented by a rather different geometry, involving a localized solitonic D3-brane. It was shown that 
the gravity solution and deconfined QCD belong to the same phase; in addition, [I] showed that the new 
geometry too offers a mechanism of chiral symmetry restoration (XSR)Q 

There have been many follow-up works on finite temperature QCD from holography which support 
the MM proposal. In particular, it is argued in j6] that the estimate of the deconfinement transition 
temperature obtained from MM fits phenomenology better than the one obtained from the Witten-Sakai- 
Sugimoto model mm • The papers mm argue that an alternative description of the high temperature 
phase of QCD is possible in a more elaborate setup using D4, NS5 and D6, D6 branes, although within the 
Sakai-Sugimoto setup the high temperature phase indeed follows MM. Besides, Ref. [T5] demonstrated the 
similarity between the confinement/deconfinement transition in the gauge theory and the phase transition 
in the MM model, and a related connection can be seen even analytically in supersymmetric Yang-Mills 
theory pjlj- More recently, Ref. |20] showed that the MM model qualitatively reproduces the instanton 

x Like in [3, here too \'SR is preceded by the deconfinement transition, although unlike in [Tj, the latter is a Gregory- 
Laflamme transition. 
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density in lattice QCD while the black D4-brane geometry does not. 

The MM construction of the deconfinement phase, however, has an unsolved question regarding the 
monodromy of quarks around the temporal circle; we will discuss this below in Section [2j The purpose 
of this article is to resolve this issue and thus firmly establish the MM method [I] for investigating finite 
temperature QCD via holography. We achieve this resolution by considering, as already hinted at in 
MM, an imaginary baryon chemical potential [21] (see Sections [3] and 0]). We find that this leads to the 
correct monodromy of the quarks suitable for describing thermal physics. As an additional application 
of our technique, we compute in Section [5] the dependence of the free energy on the imaginary chemical 
potential; we find that the result is consistent with the Roberge-Weiss potential |2~llR which is obtained 
from perturbative QCD at high temperature. 


2 Deconfinement phase in Sakai-Sugimoto model and the issue of quark 
monodromy 

To obtain S\J(N C ) QCD with Nf flavor quarks, Sakai and Sugimoto |3j considered the following brane 
configuration 


(t) 12 3 

N c D4-brane 

Nf D 8 /D 8 -brane — — — — 


(4) 5 6 7 8 9 


(1) 


Here the parentheses denote directions compactified on S 1 . To represent finite temperature, spacetime is 
regarded as Euclidean; t parameterizes the Euclidean time circle with period /3 = 1/T. x A is the Scherk- 
Schwarz circle, with antiperiodic boundary condition for fermions, which is crucial to break supersymmetry 
and to obtain the (non-supersymmetric) four dimensional QCD at low energy j2]. We define the period 
of this circle to be L 4 . 

We will consider below Nf -C N c so that we can treat the D 8 -branes as probes. In an appropriate 
large N c limit a la Maldacena mm we can replace the D4 branes by a classical geometry, and regard 
the D 8 branes as probes coupled to such a geometry. At low temperatures, the geometry, which in fact is 
the dual gravity description to compactified 5 dimensional SYM theory [23] , is given by 

3 n 


3/2 


(is 1 = 0 ! —!= ( dt 2 + V' dxj + f 4 (u)dx{] + 

.v^ V ti J u/ 


du 2 

h(u) 


+ trdP" 


hiu) = 1 - 


A5 


u' 


3/2 


1/2 


(2tt) 2 N c l ^ 5 / 4 n / 


(2) 


2 See, e.g., [22] for some recent progress on QCD with complex chemical potential. 


2 














There is also a non-trivial value of the five form potential which we do not write explicitly. Here A 5 := 
( 2n) 2 g s l s N c is the ’t Hooft coupling on the D4-brane world-volume. g s and l s = \fa! are the string 
coupling and the string length, respectively. We will also use the dimensionless coupling A ym ■= 2 A 5 / L 4 . 
The identifications t = t + /3, 24 = 24 + L 4 are implicit. 

Since the 24 -cycle shrinks to zero at u = uq, in order to avoid possible conical singularities we must 
choose L 4 as follows 


L 4 _ y/\ 5 /4ir 1/2 fo) 

2vr “ 3 0 • ( j 

With this choice, the contractible 24 -cycle, together with the radial direction u, forms a so-called ‘cigar’ 
geometry which is topologically a disc. 

Note that the above gravity solution is reliable in a regime Xym 1 in which the stringy corrections 
are suppressed. On the other hand in the regime A ym *C 1, the system flows to IR and the massive 
modes such as adjoint fermions (and scalars), together with KK modes about the 2 3 4 5 6 -circle are decoupled 


and we obtain four dimensional QCD [3] It turns out, however, that although there is no direct overlap 
of validity between four dimensional QCD at weak coupling (A ym 1) and gravity at strong coupling 
(A ym 2 > 1), it has been a highly rewarding enterprise to understand QCD dynamics by extrapolation from 
the gravity analysis (see, e.g, 0 ! for recent developments). We should note that such an extrapolation is 
analogous to the use of the strong coupling expansion in lattice gauge theory. An important consideration 
for the success of the latter expansion is that there should not be a phase boundary between the strong 
coupling theory and the continuum theory which is defined at zero lattice coupling [23]. The point of the 
MM proposal [I] was to ensure that such phase boundaries were not encountered between gravity and 
four-dimensional QCD either at low or at high temperature. 

To investigate high temperature regime using the gravity dual in the sense mentioned above, [1] 
pointed out that it was imperative to impose the periodic boundary condition along the temporal circle 


for the fermionsj 4 ! As we increase the temperature, strings which wind the temporal circle become lighter, 
and the type IIA gravity analysis is not reliable when the temperature reaches around T ~ V^ym / L 4 . 
This leads us to a type IIB supergravity description obtained by performing a T-duality on the temporal 


circle. |f| Then the analysis is reliable if T V A ym / L 4 0 . This T-duality maps the brane configuration 


3 The theory we obtain is actually four dimensional pure Yang Mills theory [5] with external quarks in the fundamental 
representation coming from the D4-D8(D8) open strings. 

4 The usual antiperiodic boundary condition leads to the black D4 brane at high temperature, which is separated from four 
dimensional deconfined QCD by a phase boundary, whereas the use of the periodic boundary condition leads to a different 
gravity solution, as we will recall presently, which does not involve such a phase boundary. The consequence of this choice 
for the description of high temperature QCD is detailed at the end of this section, and is the main issue addressed in this 
paper. 

5 If we had taken the periodicity of the fermions along the temporal circle to be anti-periodic, the system would be mapped 
to type OB rather than to type IIB. 

6 To be more precise, the mass of the winding string depends on its position, and the lightest one is located at u = uq 


3 







of the Sakai-Sugimoto model (JT|) to 


(t') 1 2 3 (4) 5 6 7 8 9 

N c D3-brane — — — — (4) 

Nj D7/D7-brane ___ _____ 


Here t' is the dual temporal circle, and its periodicity is /3' = (2ir) 2 //3 = (27r) 2 T. 

The IIB supergravity analysis shows that around 

T G l ~ ^ (5) 

-04 

the T-dual of the gravity solution ([2]) becomes unstable due to the Gregory-Laflamme instability [25], and a 
Gregory-Laflamme transition occurs. The stable solution at higher temperatures is the localized solitonic 
D3-brane geometry. Generally the analysis of the localized geometry in a compact space is difficult; 
however, an approximate geometry of the localized solitonic D3-brane geometry, can be described by the 
metric |T| 


ds 2 = a' 


H 


' 1//2 dx 2 + (1 + 2&)dx^J + H 1 ^ 2 ( 1 — ( du 2 + dt' 2 + u 2 dVL :|) 


\i =1 




2A s /p 


(u 2 + (t r — t' 0 — n/3') 2 ) 2 ’ 


= 


A5 


2irN c /3 ’ 


U TT 

$ =- 

2 


E 


m 2 + (f' — tg — n/3') 2 




— v^2A 5T 


7T 


2L 4 


( 6 ) 


The approximation used here is valid in a region far from the localized D3-branes and for T uh- Note 
that this solution describes D3-branes localized around t’ = t' 0 and it breaks the translation symmetry 
along the f'-circle spontaneously. We fix t' 0 = 0 in the following discussions. [I] argued that it is this 
localized D3-brane geometry that corresponds to the deconfinement phase in QCD, and the Gregory- 
Laflamme transition leading to this localized solution corresponds to the confinement/deconfinement 
transition^ 

To study quarks in the deconfinement phase, we need to study the probe D7-branes (jU) on the 

background geometry Q. However, we appear to encounter a problem here. As we explained above, 

we have imposed a periodic spin structure (around the temporal circle) in our supergravity analysis. 

in the IIA flame. However after the T-duality, since the radius is flipped, the heaviest one is from u = uo in the IIB flame. 
Thus, if T ^/Xym/La, while the IIA supergravity analysis around u ~ uo is not reliable, the IIB supergravity is, although 
the validity of the latter away from this region is not ensured. However, since it is only the region u ~ Uo which is involved 
in exploring the confinement/deconfinement transition, the use of IIB supergravity in T 2S> \/\ym /X 4 is justified. 

7 The spontaneous symmetry breaking of the translation symmetry corresponds to the breaking of the center symmetry 
in the deconfinement phase. 
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Naively, this would seem to suggest|j periodic BC on all QCD fermions, while the standard description of 
thermodynamics requires the BC to be antiperiodic. The thermodynamics of the pure Yang-Mills sector is 
not affected by this problem, since at weak coupling, the adjoint fermions on the D3-branes (D4-branes in 
the IIA description) are anyway decoupled from the theory, hence the choice of their temporal boundary 
condition is immaterial [T]. However, since the fundamental quarks of QCD in the Sakai-Sugimoto model 
are not decoupled in this limit, a periodic temporal BC for them would appear to be at variance with 
standard rules of QCD thermodynamics. The resolution of the above problem is a central issue of this 
paper. 


3 Imaginary baryon chemical potential and quark monodromy 

In this section, we consider the issue of monodromy of fermions in QFT in somewhat more general 
terms. Let us suppose that the quarks if) in four dimensional SU(N c ) QCD satisfy the following boundary 
condition around the temporal circlj^l, 

#* + / 3) = -e~ ie m- (7) 

As shown in Appendix [Aj a path integral with the above boundary condition can be interpreted as a 
thermal partition function with an imaginary baryon chemical potential (fi = id) 

Z{P, 9) := tr L-PH+iON^ ) ( g ) 

Here N is the baryon number operator. In Appendix [A] and what follows, we show how such a chemical 
potential arises from in the presence of a thermal Wilson line or, equivalently, from the temporal com¬ 
ponent a constant 17(1) external gauge potential. We will indicate in the next section, how such a gauge 
potential is automatically present in our problem (it is determined by the location of the D7 brane on the 
temporal circle). We will find that 

(a) Choosing the above parameter (location of the D7 brane) judiciously, we can obtain a thermal parti¬ 
tion function starting from a path integral over periodic fermions (Section [4]). This resolves the issue of 
apparent conflict between periodic fermion BC and thermodynamics. 

(b) Keeping this parameter arbitrary, the holographic setup allows us to determine the dependence of the 

free energy as a function of the imaginary chemical potential (Section [5]). The computation agrees with 

the corresponding computation from perturbative QCD— a result called the the Roberge-Weiss potential. 

8 Strictly speaking, the periodic spin structure in gravity refers to gravitinos and other bulk fermions. To connect it to 
the spin structure of the gauge theory, we need to consider the AdS/CFT dictionary which relates gauge invariant boundary 
fermion operators to bulk fermions. The naive assumption referred to here is that the spin structure of the gauge theory is 
simply the same as that of supergravity. As we will see below, in our model they differ by an imaginary baryon chemical 
potential. 

9 The temporal circle is Euclidean and is labelled by t. Note that 9 = n corresponds to the periodic boundary condition 
and 9 = 0 to the anti-periodic one. 
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4 The imaginary chemical potential in the Sakai-Sugimoto model 


We will now show how an imaginary baryon chemical potential such as in d8j) presents itself in the Sakai- 
Sugimoto model in our context (see [26] for a related discussion) @ Let us first recall that the standard 
(real) chemical potential in Euclidean SU(iV c ) QCD is equivalent to a constant pure-imaginary temporal 
component of the Euclideanised diagonal flavor gauge field. Hence an imaginary chemical potential is 
related to the real-valued Euclideanised flavor gauge field (see ([2111 for a precise relation between the 
two). In the Sakai-Sugimoto model, the flavor gauge field corresponds to a gauge field on D8/D8-branes. 
Combining this observation with the fact that fermions on the string theory side are periodic in our 
setup, we find that a thermal partition function with an imaginary chemical potential 9 can he regarded 
as a periodic functional integral with an external flavor gauge potential A{ through 


/3A{ = 9 + 7T. 


( 9 ) 


Note that the factor it comes from the difference of the periodicities. We give a detailed derivation of © 
in Appendix IA.31 

We discussed the imaginary chemical potential above (and in Appendix [A]) directly in terms of quarks 
in the D4-D8 open strings. In the dual supergravity theory, the quarks are not present directly, but 
baryons are. It is easy to demonstrate how baryons pick up nontrivial monodromies around the thermal 
circle in the presence of a constant diagonal flavor gauge field. A baryon in the gravity dual is represented 
by a ‘baryon vertex’ (a D4-brane which wraps the S' 4 of the geometry ([2])) [3j. The D4-brane couples to 
the diagonal flavor gauge field (considered as external) through a Chern-Simons term as 



The gravity functional integral in the presence of this D4-brane defines the baryon wavefunction. As one 
takes t —>• t + /3, the functional integral, and hence the baryon wavefunction, picks up an extra phase 

y 

e iN,.f,A t q ue the presence of the flavor gauge field. 

Thus, to introduce an arbitrary imaginary chemical potential, we need to simply turn on an appropriate 
constant, diagonal, flavor gauge held A {, according to ([9]). To proceed further, we need to do a bit more 
work. ([9]) is written in terms of type IIA variables. As explained in Section [2] to describe the high 
temperature limit in our holographic framework (corresponding to deconhned QCD) we need to T-dualize 
on the temporal circle. In that case, 2irA{ is mapped to the position th 0 of the D7/D7-branes on the 

10 Note that the imaginary chemical potential does not affect the periodicity of the color adjoint fermions and they remain 
periodic. 

11 In the holographic setup t'f corresponds to the point on the f'-circle reached by the D7 brane at the boundary, i.e. 
t’ f =lim u _ foo t'(u). 
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(11) 


r@ 


dual temporal circle. The type IIB version of (J9l) is, the 

t 'f = J^ + 6) = ^ + 0). 

In particular, we can satisfy the above condition with (this was indeed the choice made in [1]) 

t'f = P'/2 (=» 0 = 0). (12) 

Recall that by convention we have put the D3 branes at t' = 0. Equation (11211 implies that if the D7/D7- 
branes are placed antipodally from the D3 branes on the temporal circle of the IIB theory (equivalent to 
a choice 0 = 0), the path integral with periodic fermions does reduce to a thermal partition function, thus 
resolving the apparent puzzle raised in Section [2j 

One might wonder whether we can place the D7 branes at will at any point on the temporal circle. In 
case the D7 branes are infinitely heavy, the position of the D7’s can be regarded as an external parameter. 
However, in the next section, we will consider the potential energy of the D7 branes, placed at a position 


t'f at the boundary (u oo) ^ 


of the potential 


0 


We will find that the choice 0 = 0 mentioned above is indeed a minimum 


5 Free energy as a function of imaginary chemical potential 


In this section we will calculate the free energy of the D3-D7 system (primarily as a function of the 
location t'f of the D7 branes on the temporal circle). As mentioned before, we will use the dual gravity 
description Q for the D3 branes and consider the D7 branes as probes coupled to this geometry through 
a DBI action. In the deconfinement geometry ([6|), there are two possible stable configurations of the 
D7-branes which correspond to the chiral symmetry being preserved/broken [lj (see Fig. [[]). 



Figure 1: D7/D7-branes in the localized solitonic D3-brane background ([6]) . There are two possible 
configurations of the D7/D7-branes: (a) the D7 and D7-branes extend parallel to each other, which 
implies that the chiral symmetry is preserved, (b) the D7 and D7-branes are connected, which implies 
that the chiral symmetry is broken. 


12 The length of the temporal circle in the IIB is denoted by /3' which equals 4tv 2 //3. 

13 This can be regarded as the potential energy between the D7 and D3 branes in the open string description. 

14 See (1181) . Other minima, corresponding to non-zero values of the integer k, correspond to the presence of k D3-brane 
pairs as explained in footnote [25] 
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Let us consider the chiral-symmetric configuration for simplicity; it is straightforward to generalize our 
considerations to the broken symmetry configuration. We assume that the branes are not curved in X 4 
space. Then the induced metric on the D7-branes is given by 


^ S D7 =a ' 


H ~ 1/2 ( dxj ) + LT 1/2 


v 1=1 


1 +[^M) ] du 2 + u 2 dQ 2 d 

du 


n 

_ A 5 
47r 


2A5//3 


(u 2 + — n/3 ') 2 ) 2 


Re 


1 


tanh 


7 r(u—it'(u)) 

V 


+ WT2 Re 


1 


f3'u 


2 / n{u—it , (u)) 


sinh 2 


P' 


(13) 


where t'(u ) is the position of the D7-brane subject to the boundary condition t'(u —>• 00 ) = t’j. We have 
assumed that the D7-brane is far away from the D3-brane and neglected the quantity 4? appearing in the 
metric ©. Later we will justify the validity of this assumption. The DBI action for the Nf D7-branes 
then becomes 


cD7 

°DBI 


= N f T 7 I d 


/■ 


xc ^ < ^°°) 


\J det 


9d 7 


8 vr 2 N f N c B r 

3(2vr)6A 5 3 Jo 


duu^H 1 ' 2 




(14) 


Here V 3 is the volume of the xi,X 2 ,X 3 space, and we have used e^°° = g s and T 7 = 1/g s (2ir) 7 a ' 4 which 
is the tension of the D7-brane. The Euclidean action has a minimum at the constant value t'[u) = ft'/2 
(mod /?')@ Thus the stable configuration of the D7-branes is t'{u) = t'j = j3'/2 with dt'{u)/du = 0. This 
configuration is precisely (11211 which corresponds to 9 = 0 @ 

To calculate the t'j dependence of the DBI action, we approximate dt'(u)/du = (j^and take t'(u) = t/j. 
The integral (fTTl) evaluates as 

J duu 4 H 1/2 (c 0 + c 2 (Pt' f - 2tt 2 ) 2 -i - ), (15) 

where cq is a divergent constant which should be removed through a regularization and C 2 = 20.7022... . 

This gives us the potential of a single D7 brane as a function of the position variable t'^. As we found in 

Section U is directly related to the imaginary chemical potential 6 , according to the relation dill) . By 

substituting this relation into the above and summing up the contributions of all D7 and D7-branes, we 

15 To prove this, note that the minimization of the action requires (a) setting dt'/du = 0 and (b) minimizing H. From the 
second line of m H, for constant t'(u) = t/ is clearly minimized at t/ = /3' /2 (mod /3'). 

16 Ref. ]T] investigated the chiral symmetry breaking/restoration transition by comparing the DBI action for the configu¬ 
rations of the D7-branes corresponding to these phases (see figure [T]). In this study, the D7-branes are set at t' = /3'/2, since 
this is stable for both chiral symmetry preserved and breaking configurations. Through the relation CD, it corresponds to 
9 = 0. This shows that the results in [T] are valid for the standard thermal quarks. 

17 If we do not use the approximation dt'(u)/du = 0, C 2 is modified but qualitative properties would not change. 


















obtain the 0 dependent classical action 


C-D7D7 

'-’DBI 


2N c N f V 3 n2 , 

SvaC ® + "" 


(16) 


Once again, we see that the choice 0 = 0, made in pQ, and given by () 12 H is a consistent classical solution 
according to this potential. 


This, however, is not the end of the story. The partition function of QCD with an imaginary chemical 
potential © has a discrete symmetry 


Z{T-6) = Z{T-9 + 2n/N c ), 


(17) 


which is a generalization of the Zyv c symmetry of pure Yang-Mills theory to QCD (see Appendix I A. 41) . 
The potential (fT6l) clearly does not respect this symmetry. To fix this problem, we need to understand 
how the above Zyv c symmetry is realized in the gravity dual. We find, following similar considerations in 
E3 E3 ESI ESI GEI , that the missing pieces of the puzzle are provided by the dynamics of some bulk gauge 
fields, namely, in our setup, the graviphoton (KK gauge held) and C 4 gauge held. We give the details in 
Appendix [Bj the upshot is that the dynamics of these gauge helds changes the classical action (11611 into 


S(9) 


2 N c N f V 3 ^ 

3vr 2 \/WLZ 


C 2 min 
fee z 



(18) 


The above action, plotted in figure El obviously has the discrete symmetry 9 —> 9 + 2ir/N c . 



Figure 2: A schematic plot of S(9), the effective potential for the imaginary chemical potential 9. 


The potential has N c minima for 9 (0 < 9 < 2ir) and indicates N c phase transitions as 9 varies. Note 
that the minimum discussed above, namely 0 = 0 (corresponding to (|12p f remains a minimum for the 
above potential as well. Through each transition, the distance between the D7-branes and the D3-branes 
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is shifted by /3'/N c [J£| and they are always antipodal at the minima of the potential. It ensures the 
assumption below (fT3l) . 

In [21] i Roberge and Weiss obtained the following result for the potential for an imaginary chemical 
potential 8 in the deconfinement phase from perturbative QCD: 


Wert. {8) = min 
p fcez 


12 M vr 2 \ N c 


(19) 


Remarkably, this potential qualitatively agrees with the holographic result (1181) we found above. 
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A A review of imaginary chemical potentials in field theory 

We summarize properties of the imaginary chemical potential which are used in the main part. 

A.l Relation between imaginary chemical potential and periodicity of fermions 

First we discuss the connection between the imaginary chemical potential and the periodicity of the 
fermions. In general, for a system described by a Hamiltonian H with fermions described by dynamical 
variables rfi (where the index i refers to any set of indices of the dynamical variables including the 
coordinates), we can derive the following identity 


tr ( e -Aff+i0iV) = f 

Jfj 

-I 


[dfj*dfj\ exp < — f dt 

v(8)=-fK o) I Jo 


E 


r)*d t r\i + H(ff*,ff) - i/3 




[dfj*dff\ exp { ~ I dt 


l 


^2rn*d t r]i + H(ff*,ff) 


. . ( 20 ) 

/ rf(/3)=— e~ l6 fj{0) 

where t is the Euclidean time and N := Y2iVi*Vi is the fermion number operator. The trace on the left 
hand side is taken over the fermions. We can derive this relation in a usual manner with the slicing in 
/3-direction and the coherent state representation. The nontrivial periodicity on the second line results 
from the identity e %eN \ff) = \e l6 ff) on the coherent state | if) in the trace. Therefore we can identify the 
imaginary chemical potential with the monodromy of the fermions around the Euclidean time circle] 19 


18 To be more precise, through each transition, tj in dm is shifted by 1/V C , and hence t'f is shifted by /3'/N c due to the 
relation above (S3. 

19 Two known special cases of l|20|l are worth noting: (a) a thermal partition with no imaginary chemical potential (6 = 0) 
corresponds to an antiperiodic boundary condition around the Euclidean time circle, and (b) a thermal partition function 
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A.2 Relation between imaginary chemical potential and gauge fields 

Next we establish the connection between the imaginary chemical potential and the U(l) component of an 
external gauge field. Let us consider a QCD-like theory with a color gauge group G c and matter fermions 
ip in the fundamental representation of G c and some representation of a flavor group Gf. We gauge both 
G c and Gf with the gauge potentials A c and respectively, and regard A c as dynamical and Af as 
external. The Lagrangian density (in the Lorentzian signature) is 

= -ip^D^ + m)ip 


with the following definitions: 

D^'-d^-iA^-iAfy, {'f, r f} = 2rf v , rj^ = diag(-l, 1,1,1). 

The Hamiltonian reads 

H m [A f ] = J d 3 x [~ip^A 0 ip + + rmpip], 

where A := A c + A*, and we displayed the external flavor field A? explicitly in the Hamiltonian. Note 
that the index 0 stands for the Lorentzian time while t stands for the Euclidean time. 

Now let us turn on only the temporal component Aq of the external flavor gauge field A? and suppose 
that this Aq is proportional to the identity matrix. Then the first equality of (1201) with the Hamiltonian 
H m [Af] and the number operator N := f d 3 xip^ip reads 


tr exp 


I ip(f3)=-'ip(0) 


-PH[A f 0 \ + i6N 

[dA c d^* di’} e~ s ^ AC 


exp 


rP _ 

— / dtd 3 x [ip^dtip + ip(fy l Di + m)ip — iift A t ip — i / 3 _1 dip^ip] 
Jo 


f [dA c dip*dip] e exp | — f dtd 3 x [ip^dtip + ip(ff r Di + m)ip — iip^(3 1 (/3A t + 6) ip] 
/^(/3)=—-0(0) l Jo 

-PH[0 \ + i(6 + PA{)N 


= tr I exp 


( 21 ) 


where we defined At := — iAo=A c t + A^t, and «S , C [^4 C ] is the Euclideanized Yang-Mills action of the color 
gauge field The final expression makes it clear that A{ (times j3) plays the same role of an imaginary 
chemical potential 9. 


with an imaginary chemical potential 9 = n corresponds to the periodic boundary condition. The latter observation ties up 
with the fact that an insertion of e l6N in the LHS of (1201) with 6 = 7r is equivalent to inserting (—1)^, which, of course, 
converts a thermal partition function into a Witten index, and is hence represented by a path integral over periodic fermions. 
20 The trace is taken not only over the matter fermions but also over the color gauge field. 
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A.3 Derivation of OH} 


We will find a string-theoretical setup which yields the thermal partition function with an imaginary 
chemical potential. The thermal partition function is given by a path integral of a Euclideanized QCD- 
like theory with anti-periodic fermions while a gauge theory on the D4/D8/D8 branes in our string- 
theoretical setup is a Euclidean QCD-like theory with periodic fermions. If we add a diagonal gauge field 
At on D8/D8, this couples with the fermions as a flavor gauge field. Thus the string theory side gives us 
the following path integral^! 


[ [dA c dip*dip] e 5c ^exp|— f dtd 3 x [ip^dtip + ip(rfDi + m)ip — iifij3 ■ 

J^P)=m {Jo J 

By the second equality of PUD , this path integral is equal to 

f [dA c dip*dip\ e~ Sc ^ A J exp/— f dtd 3 x [ip^dtip + ip{^ l Di + m)ip — iip^ /3~ 1 (f3At — 7r)?/>]l 
J V’(/3)=—b(o) {Jo J 


( 22 ) 


(23) 


This is equal to the thermal partition function with an imaginary chemical potential /3At — it. In other 
words, if we want to compute a thermal partition function with an imaginary chemical potential 6 , then 
we may consider a D4/D8/D8-brane setup with a diagonal gauge field /3~ l (9 + n) on the D8/D8-branes. 


A.4 Periodicity of thermal partition function with imaginary chemical potential 

Finally, we explain the periodicity (1171) of the partition function Z(T,6) := tr{e^ l3H+teN ), where G c = 
SU(A C ) and Gf = U (Nf). Its path integral expression can be written down using (l2T|) as 

Z(T,9)= f [dA c di/j*dip] e- s ^^’ AC] - Sc[FC \ (24) 

rP _ 

S[ip,ip,A c ] := / dtd 3 x [ip^dtip + ip^Di + m)ip — i^A^ip]. 

Jo 

Let us consider the following SU(A C ) transformation g with a Zjv c -twisted boundary condition 

5 (/3, x) = e - 2 - fc /%((), *), g(t, x) G SU(A C ), k = 0,1, ■ ■ ■ , N c - 1. (25) 

Under this transformation, the fermions ip and the Polyakov loop W ( A c ) transform as 

ip(t, x) ^ ip°(t, x) = g(t, x)ip{t, x), W(A C ) i y W(A C9 ) = e~ 2irik / Nc W(A c ). (26) 

This transformation changes the periodicity of the matter fields as follows 

4’ 9 (P,x) = -e~ ie iP( 0,*) —> ip 9 (P,x) = -e-^ +2nk / Nc) iP 9 {0,x). (27) 

21 More precisely, a classical supergravity plus the D8/D8-brane actions evaluated at an on-shell configuration gives us 
(12211 of the corresponding gauge theory, which is a QCD-like theory with periodic fermions. We perform this in the T-dual 
setup in Section [5] 
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Then we can derive the periodicity (1171) as follows: 


Z(T,9)= I 

J %j) 

-f 

J 'll) 


i/j(/3)=— e -lS i/j(0) 


[dA c dil)*dif>] e -s[M,A]-Se[F*] 


^9^) = - e -i(S+2-Kk/N c )^g^ 

= Z(T,d + 2irk/N c ), 


[dA^d^dr] e sir^,A^]-S C [F^] 


(28) 


where in the second equality we used the invariance of the path integral measure and the action under g. 


B Ztv c symmetry for 0 in gravity 

In this section, we show how the discrete Zn c symmetry 6 —>• 9 + 2n/N c discussed around (1171) is realized 
in the localized D3-brane system. The origin of this symmetry, as noted above, is the fact that the color 
gauge group dual to the near horizon geometry is SU(N c ) (rather than U(N c )) with a center Z?v c . Related 
symmetries in black brane geometries have been discussed in [271 [30] , and our derivation is similar to these 
works. 

To see the Zjv c symmetry in type IIB supergravity, we consider a Kaluza-Klein decomposition on 
the f ; -cycle and study the low energy dynamics in terms of the 9 dimensional IIB supergravity. The 9d 
supergravity action has the following termini 


P 


17^/4 \ 2 


- / d y x y/g F£ + / Ai A F A A F 4 , 


(29) 


(27r) 7 Q: , ' 

where T 4 = dC%, the three-form potential C 3 is the dimensional reduction of the four-form potential C 4 , 
and the one-form potential A\ is the KK gauge held from the 10 dimensional metric. We will show that 


r) 


1 

= pJu 0 


A u du 


(30) 


takes a discrete value k/N c for some integer k at low energy— a fact which, we will see, reflects the Zjv c 


symmetry]^ 


At low energies, the kinetic term of (|29l) reduces to 


2e 2 Jr3xS} 


d 4 x Fj 


(31) 


22 These are related, by T-duality along the d-cycle, to the terms 

| J d 10 x^g Ft + j B 2 A F 4 A F 4 

in IIA supergravity which are employed in EDI to explore the Zjv c symmetry in the black D4-brane background. 

23 mo denotes the IR end of n-coordinate in the 9 dimensional theory which is related to the position of the horizon of the 
localized D3-brane geometry. Although the position of the horizon depends on t' and u in the original 10 dimensional system, 
within the KK-decomposition the t' dependence of the metric is mapped to the condensation of KK non-zero modes. 
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through the integration over S ' 4 and the u direction. Here e 2 is an effective 
(1291) becomes 



The CS term of 


F 

(2n ) 7 a’ A 


Hi A F 4 A F 4 



(32) 


where we have used fg4 F4 = N c and fj ,3 := (27r) ‘W 2 is the tension of a D3-brane. 

the Hamiltonian density from the action (1291) reads 

u = y (n - fi 3 p' V N c ) 2 , 


Therefore 


(33) 


where n := F/^/e 1 + fi^r]N c = —id/dCs is the conjugate momentum of C 3 with regard to one of the three 
directions (xi,X 2 ,X 3 ) of R 3 (see O) as the ‘time’ direction for the canonical quantization. Then similarly 
to 2d QED [3D [32], a basis of energy eigenstates of this system can be taken as [23130] 


with energy eigenvalue 





e 2 /i|/T 2 


(k — rjN c ) 2 


(34) 


(35) 


Here k must be integer-quantized. Note that under a shift r/ —>• r]+l/N c , there is a monodromy Tfc —» 
accompanying a pair creation of D3-branea 25 l (this is an analogue of a pair creation of the electrons in 2d 
QED [31]). Therefore the potential for rj becomes 


eV !/? /2 


V ( 77 ) = min 
fcez 


(k - r]N c 


(37) 


and is minimized at rj = k/N c for some integer k. Hence ij is restricted to discrete values at low energy. 

Now let us show that the shift tj — > r]+l/N c induces the Ztv c symmetry for 9. We have considered the 9 

dimensional supergravity so far and ignored the KK non-zero modes. However they should be included in 

24 Although u integral is over an infinite range, the effective coupling e 2 is finite, as in the case of a similar discussion for 
black branes }27i |30| . 

25 The brane configuration of the created D3-branes are given by 

(t') 1 2 3 (4) 5 6 7 8 9 

N c D3-brane (LSD3) — — — — (36) 

pairs of D3-branes — — — — 

where we have taken x 3 as the ‘time’ direction for the canonical quantization. Note that this brane configuration is the 
T-dual of the black D4-brane system studied in [30l . 
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the system, and importantly they always couple to the KK gauge field A u through the covariant derivative 


d u i 


27771 




(38) 


where denotes the n-th KK mode of field (j>( = e 2mnt ^ <f> n ) which symbolically denotes all fields in 
the nine dimensional supergravity except A u . Now A u can be absorbed into the phases of the KK modes 
by the field redefinition 


= exp 


. 2 irn 

! y 


du' A u ) </>. 


1 U o 


(39) 


Using this fact, we show that we can read off the solution of supergravity with non-zero A u from the 
solution with A u = 0. Since A u with no derivatives appears only at the CS term (1291) after the field 
redefinition (1391) , the equation of motion for cj) n is the same as one in A u = 0 case except for the RR fields 
which couple to A u directly through the CS coupling. Thus if we ignore the backreaction of these RR 
fields, 4> n = 4 > n \ A =Q is a solution of supergravity with A u / 0, where </> n U =o a solution of supergravity 
with A u = 0. Therefore in terms of the original ten dimensional coordinates 


,, , , V- (2TTnt' . 27T71 f u . , . \ , . 

^,a) = 2 J ex P^—+iy j du A u j (/> n \ Au =o- 


(40) 


is a solution of the supergravity with non-zero A u . It means that the solution is obtained just by replacing 
t' ~>t' + f Ug du'A u . 

Let us consider how the free energy for the imaginary chemical potential (1161) is modified by A u . The 
free energy depends on 9 which is related to the position of the D7-branes at the boundary {u —> oo) 
via (HU). Thus according to the above arguments, we replace t f —> t' + duA u = t! + f3'ri (and hence 
9 — >• 6 + 2 irr] by (fill) ) to obtain 


SBS ? 7 = 2 -MAL rh 2 (9 + 2n v f . 


2>TT 2 y /Ay M L 4 

By adding the energy (|37l) for F 4 , we obtain the ^-dependent terms of the classical action 


SR™ 7 + S F , = tIc 2 (9 + 2 + mill (k _ vNc f 


-’DBI 


(41) 


(42) 


37T 2 y/ XyM L 4 ' ' fcGZ 2 

Then after an integral over rj the system acquires the symmetry 9 —> 9 + 2tt/N c by compensating the shift 
of rj as expected in QCD (1171) . By regarding this symmetry, we can rewrite the 0-dependent part of the 
classical action 


S(9) = 


2N c N f V 3 ^5 


3ir 2 \/ XymLi 


T 2 C 2 min 9 — 


fee z 


2 -k k 


(43) 


Thus we obtain (fT8l) . 


26 We performed an integral over r /, observing that the coefficient of 77 2 in the second term of (1421) is much larger by N c 
than that in the first term. 
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